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Abstract 

In the effective theory for the deconfining phase of SU(2) Yang-Mills ther- 
modynamics we compute estimates for the moduli of the irreducible three-loop 
diagrams contributing to the pressure. Our numerical results are in agreement 
with general expectations. 



1 Introduction 



To obtain essential analytical insights into the thermodynamics of four- dimensional 
Yang- Mills theories is difficult even at high temperature where gauge fields do propa- 
gate. A perturbative treatment of Yang-Mills thermodynamics, which is technically 
highly involved, runs into problems that are associated with the masslessness of the 
fundamental, propagating degrees of freedom. A manifestation of this problem is the 
apparent nonconvergence of the perturbative serie^ which can be pinned down to 
the weak screening in the magnetic sector of the theory [3]. Loosely speaking, this 
problem arises because the perturbative a priori estimate for the thermal ground 
state is inappropriate. Namely, the presence of topologically nontrivial fiuctuations, 
which do contribute to the thermodynamics of the Yang-Mills systemic in a direct 
(ground state) and indirect (masses) way, is neglected in perturbative loop expan- 
sions. 

In [6] a nonperturbative approach to SU(2) and SU(3) Yang-Mills thermody- 
namics is developed. Let us discuss the deconfining phase for the SU(2) case only. 
The idea is to first derive a thermal ground state which is composed of interacting 
topological fluctuations: calorons and anticalorons. While the dynamics of these 
fluctuations would appear to be highly complex in the hypothetic case, where an 
externally provided probe of a given momentum transfer is applied to the system, 
the case of pure, inflnite-volume thermodynamic^ selfconsistently adjusts a maxi- 
mal resolution |0| in such a way that the ground state admits a remarkably simple 
analytical description. To derive this situation a spatial coarse-graining needs to be 
performed which yet turns out to be sufficiently local to only require consideration 
of calorons and anticalorons of topological charge modulus |(5| = 1 [7]. As a conse- 
quence, an adjoint and inert (spatially homogeneous) scalar field (p emerges which 
together with a pure-gauge (coarse-grained) configuration describes the ground state 
of the system, see also [S]. The presence of signals a dynamical gauge-symmetry 
breaking SU(2)^U(1) implying that two out of the three propagating and coarse- 
grained gauge-mode species acquire a (temperature-dependent) mass. This and the 
fact that the off-shellness of these modes and the momentum transfer in local ver- 
tices is highly constrained after spatial coarse-graining imply the rapid convergence 
of loop expansions in the effective theory [9J. The purpose of the present article 
is to demonstrate this by estimating the irreducible three-loop contribution to the 
pressure in deconfining SU(2) Yang-Mills thermodynamics. 

The paper is organized as follows. In Sec. [2] we briefiy review technical essentials 
of the effective theory. Sec. [3] first explains what is meant by the term 'irreducible 

"'^At T = it is suggestive that the perturbative series represents an asymptotic expansion, and 
so the first few orders do capture a great deal of the physics [U |2] • The perturbative series at finite 
temperature apparently does not enjoy this property. 

^For example, the trace of the energy-momentum tensor rises linearly with temperature [3], and 
the dimensionally reduced theory confines [5]- 

■^Subject to two scales only: temperature T and Yang-Mills scale A. 
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diagram' according to the discussion in [9j . Subsequently, we elucidate the structure 
of irreducible three-loop diagrams. In Sec. H] we perform the integrations numerically 
using the Monte-Carlo method. Results are indicated and discussed. Finally, SecO 
gives our conclusions. 

2 Effective theory for deconfining phase 

Here we very briefly review the effective theory for SU(2) Yang-Mills thermodynam- 
ics being in its deconfining phase. The following effective action emerges upon a self- 
consistent spatial coarse-graining involving interacting nontrivial-holonomy calorons 
of topological charge modulus |Q| = 1 (ground state) and topologically trivial gauge 
fields (excitations) [6]: 

S = tT j^^rj Q G^,G^, + D^(t>D^(t> + A ■ (1) 

In Eq. (HD G^, = G'^.f , = d^A'i,-d,A- + ee-''^AlAl, and D^<p = d,<P + te[<P, A^] 
where A^ is the (coarse-grained) gauge field of trivial topology, and e denotes the 
effective gauge coupling. 

In a first step, only a noninteracting (trivial-holonomy) caloron and anticaloron 
is coarse-grained over an infinite spatial volume into the phase of a spatially homo- 
geneous adjoint scalar field 0. Subsequently, the presence of a Yang-Mills scale A is 

assumecElto obtain the modulus \<p\ = \J The latter determines the finite spatial 

length scale at which the above coarse-graining saturates. Taking interactions 
between calorons and anticalorons into account, a (coarse-grained) pure-gauge con- 
figuration emerges. Going from a gauge, where winds along the compactified 
euclidean time dimension, to unitary gaug^ we arrive at 

>Cr- = ^ M = \ {GTK]f + 2e' 101^ {{alY + [alf) + 2^ . (2) 

The effective gauge coupling e enters into both the effective field strength G"^^" and 
the mass m for the fields a];"^ . One has 

= m{Tf = ml = ml = 4e^ \(f)f , rug = . (3) 

The gauge mode aj^ stays massless on tree-level in the effective theory (adjoint Higgs 
mechanism). The associated unbroken U(l) gauge freedom is fixed by imposing the 
Coulomb condition dia^ = 0. This corresponds to a completely fixed and physical 
gauge (unitary-Coulomb gauge). 

^This assumption is actually redundant since A can be interpreted as a nonperturbative inte- 
gration constant, see [5]. 

^ This involves an admissible electric center transformation [6l [10] . 
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Let us now give the one- loop expressions for the partial pressures Pi = P2, and 
P3 as exerted by the fluctuating gauge modes aj^yoj^, and aj^, respectively. One has 

90 ' io 27r2 V / 

Notice that in deriving Eq. (jl]) the vacuum part in the one-loop expressions for the 
partial pressures can safely be neglected by virtue of momentum constraints in the 
effective theory [6], see also below. The total one-loop pressure Pi.ioop (exerted by 
fluctuating modes) is given by Pi.ioop = Pi + P2 + -Ps- 

The temperature evolution of the coupling e is determined by the relation 

a = 27reA-3/2 (5) 
and by the (inverted) solution of the (one-loop) evolution equation 

where 

D{a)= / dx , ^ , (7) 

Jo vx^ + a? exp(v2;2 -\- a?) — 1 

A = and a = Eq. ([6]) guarantees the invariance of the Legendre transforma- 
tions between thermodynamical quantities when going from the fundamental to the 
effective theory. The evolution of e with temperature exhibits a logarithmic pole, 
e oc — log(A — Ac), where Ac = 13.89 denotes the critical value of the (dimensionless) 
temperature, and a plateau setting in for A slightly larger than Ac. The value of e 
at the plateau is e = -^/Stt ~ 8.89. 

In calculating radiative corrections to the free-quasiparticle (one-loop) pressure 
in a real-time formulation the loop momenta in the effective theory are subject 
to constraints. The latter emerge due to the existence of a maximal scale of 
resolutioii^. In [9] these constraints were discussed, and we content ourselves with 
simply quoting them here. First, any propagating gauge mode with four-momentum 
p cannot be further off its mass-shell then That is 

\p^ — m^\ < 101^ (for a massive mode) , \p^\ < |0p (for a massless mode) . (8) 

Second, the momentum transfer within a vertex needs to be constrained. For a three- 
vertex (ii) is already contained in (i) by momentum conservation in the vertex. For 
a four-vertex one needs to distinguish s, t, and u channels in the scattering process. 



^ Since the length scale \4>\~^ is deep inside the saturation regime for the spatial coarse-graining 
[3I2] physical quantities should not depend on a mild rescaling of |0|. We have checked the validity 
of this assertion when computing the polarization tensor of the massless mode |10| . 
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Suppose that the ingoing (outgoing) momenta are labeled by pi and p2 {ps and 
Pi = Pi + P2 ~ Ps) ■ Then the following three conditions emerge 

\iPi+P2y\ < (s channel) |(P3 - Pi)^| < |0p , (t channel) 

|(P2-P3)'| < 101', channel). (9) 

For a three-vertex conditions (Q are already contained in ([8]) by momentum conser- 
vation in the vertex. Notice that the three conditions in Eq. (Q reduce to the first 
condition if one computes the one-loop tadpole contribution to the polarization ten- 
sor or the two-loop contribution to a thermodynamical quantity, say the pressure, 
arising from a four- vertex [10] . Namely, the t-channel condition is then trivially 
satisfied while the w-channel condition reduces to the s-channel condition by letting 
the loop momentum k — > —k in — A;)^| < see [HI [TOl [H]. Notice also that 
upon a euclidean rotation po ipo the first condition in(IH]) goes over in 

\p^ + m^\ < . (10) 

For SU(2) the quasiparticle mass is given as m = 2e with e > ^/Sn ~ 8.89 [H]. 
Thus condition (fTOl) is never satisfied, and massive modes propagate on-shell only. 
Conditions ([8]) and iQ imply that the higher the loop order the more suppressed their 
contribution to a thermodynamical quantity. General arguments suggest that, apart 
from diagrams associated with one-particle irreducible resummations of propagators, 
only a finite number of diagrams contributes to the loop expansion pj. The main 
purpose of the present work is to demonstrate the validity of this on the three-loop 
level. 



3 Irreducible three-loop diagrams 

We are interested in an estimate for the modulus of each irreducible three-loop 
diagram contributing to the pressure. By three-loop irreducible we mean that the 
diagram does not include any line that is dressed by (multiple) insertions of one-loop 
polarizations. These one-particle reducible contributions to the propagator must be 
resummed to avoid the occurrence of pinch singularities [HI [lO] . Apart from a mild 
modification of the tree-level propagator this modifies the dispersion law of the 
associated mode which in turn leads to a slight modifications of the constraints in 
dH]) . The claim of [9] is that the loop expansion terminates with respect to irreducible 
diagrams, that is, with respect to all those diagrams which do not yield a one-particle 
reducible diagram upon performing a cut (in all possible ways) on a single line. 

The only three-loop irreducible diagrams are depicted in FiglH In the following 
we use the convention as in Fig[T] for labelling the loop momenta. For diagrams 
A, B, and C the number K of independent, potentially noncompact loop variables 
{Po, IpDj! (^ = 1)2,3), is = 6, and the number K of independent constraints is 
K = 7. This implies that the support for the loop integrations is either compact or 
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P4= P1+P2-P3 



P4= P1+P2-P3 
P5=Pi-P3 
B 



P4= P1+P2-P3 



A 



C 



Figure 1: Irreducible three-loop contributions to the pressure. Solid (dashed) lines 
are associated with the propagators of massive (massless) modes. 

empty, see also the discussion in [9J. As we shall see, the former possibility applies 
to diagrams A and B while diagram C vanishes. 

Let us first discuss diagrams A and B. Upon use of the Feynman rules, see [10], 
considering the symmetry factor ^, by appealing to the triangle inequahty, using 
the fact that the modes a^'^ propagate (on-shell) thermally only, integrating over 
the time components of the independent loop momenta (momentum conservation), 
and after a rescaling of the radial components of loop momenta as 



one arrives at the following estimate for the moduli of the pressure corrections 



I - Pi ( 

Pil Xi = — , ( 



1,2,3) 



(11) 
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APa,APb 



\APa{b)\ < ^ [ dxi [ dx2 [ dxs [ dzi2 [ dz^ j dz23 

6XZ X [Z7l) ^^^^^^^J J J J J Jz,,,, 

1™2 ^2 ^2 
12 3 



v/(l - zf,){l - zf,) - (Z23 - 212213)2 ^xl + 4eVa^l + 4eVa^i + 
5 (^4e2 + {-iy+"'^xl + Ae^^xl + Ae^ - X1X2Z12- 

((-l)'+ya;2 + 46^x2 + 4e2 - xiX3^i3)- 

+ 46^X1 + 4e2 - X2X3^23)) X 

(^27^X-^/^^Jxl + Ae^^ x 



VA{B){^,2,l,m,n) 



n.B \2'k\-^/^\IxI + 4e2 ) ( 2u\-''^\l xl + Ae^ \ x 



+ 4e2 + (-1)"^ + 4e2 + (-1)" Jx^ + 4e 



(12) 



where ^12 = cosZ(xi,X2), ^13 = cosZ(xi,X3), and 223 = cosZ(x2,X3). The 
functions Va-, Vb emerge from Lorentz and color contractions and are regular at 
Xi = X2 = X'i = Q (mass gap for a^'^). We refrain from quoting them here explicitly 



for environmental reasonqj- In addition, we define: 

2;23,n = COS |arccos2;i2 — arccos2;i3| , Z2-3,i = cos |arccos^i2 + arccos^isl . (13) 

The integrations in Eq. (|T2l) are subject to the following constraints (see (IHl) and 
(ED): 

Z12 < ^ (^Ae^ - ^xl + Ae^^xl + Ae^ + = duixi, X2) , 

213 > ^ (^-4e2 + .^xl + Ae'^^Jxl + Ae'^ " ^) ^ ^^^(a^i, X3) , 

223 > ^ f -4e2 + Jxl + Ae'^Jxl + Ae-' - = to(x2, X3) . (14) 

3^23^3 \ V V ZJ 

Notice that for diagram B the constraint (IE]) for the momentum P5 = pi — ps of the 
massless mode is the same as the t-channel constraint for the four-vertex. Therefore 



^Upon request the reader will be provided with the Mathematica notebooks containing the 
functions Va, "Pb in explicit form. 
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no extra condition for the (off-shell) momentum is needed. According to the 
investigation in ^ the conditions (IT^ together with Eq. (IT^ imply that the support 
for the integration in Xi,X2, and ^3 is contained in the compact set {xi, X2,X3 < 3}. 

Let us now turn to diagram C. We first consider the case that the momenta pi and 
P2 of the massless modes, compare Fig.[T], are both off-shell within the constraints 
dictated by ([8]). In analogy to diagrams A and B one then derives that 



\AP, 



c\ 



< 



3 X 25 X (27r)8 ^ 

^ ' Lm=l 



/ dyi / dxi / dx2 / dx^ / dzi2 



dz^-i I dz- 



™2 /yi2 

12 3 



- ^?2)(l - Zl^) - (2^23 - ^122:13)^ 

TiB l'27rA-3/2 

( 27rA-3/2 /x2 + 4e2 " 



Vcix,z,yi,l,m)\ X 



23 



223, i 



-l)V4 + 4e2 + (-l)™/2(x,z) 



/2(x,z)^x2 + 4e2 



(15) 



where 



/2(x, z) = Y a;^ + 2x1X22:12 - 2x1X3^13 - 2x2X3^23 , 



(16) 



and ^23,/! 2;23,„ are defined as in Eq. (fT3l) . The function Vc emerges from Lorentz 
and color contractions and is regular at Xi = X2 = X3 = (mass gap for aj^"^). The 
integrations in Eq. (fTSll are subject to the following constraints 



1 > + ^2 - a^i - 3:2 + 2yij/2 - 2x1X2212! , 



1 > |?/2'- 2^2 + 4e2-(-l)'22/2Vxi + 4e2 + 2x2X32:23 



1 > |?/?-x? + 4e2- (-l)'22/n/x2 + 4e2 + 2x1X3^13 



1 > \yt 



X 



11 ' 



1 > I1/2 



X 



21 ) 



(17) 



where 



y2 



-yi + 2(-l) 



x: 



+ 4e2 + (-l)™/2(x,z) 



As we shall see in Sec. HI the constraints in (fTTl) imply that the support for the 
integration in Eq. (fT5|) is empty. As a consequence, the cases that one or both of 
the massless modes in diagram C propagate on shell also have an empty support. 
This is because the conditions \p1\, \p2\ < which went into (fT7|) . contain the 
cases Pi = and/or P2 = 0. Thus diagram C is the first example of an vanishing 
irreducible diagram in the loop expansion. As was argued in |9] one expects that 
the number of such cases will drastically increase with increasing loop numbers. 
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Figure 2: An upper estimate for the modulus of the pressure contribution jAP^I 
due to diagram A in Fig.[Tl The plot shows this estimate normalized to the one-loop 
result Pi_ioop, see Eq. (j4]). 



4 Numerical evaluation and results 

Here we present our results obtained by using the Monte-Carlo method of integration 
for the regular integrand^ of Eqs. ( fT2i) and ( |T5l) . Let us first discuss diagrams A 
and B. It is known [_9j that the support in X2,X3 for the integral in Eqs. ( |T2l) is 
contained in the compact set {x2,X3 < 3} while the support for the integration in 
Zi2, -213, naturally is contained in the set {-1 < 212, 2:13 < +1; -223,/ < 2^23 < -Z23,n}, 
see Eq. f ll3p . Points are thus chosen randomly in the union of these two compact 
sets. Any point that satisfies the constraints (IT^ contributes to the integrals. We 
have worked with a sample size of 5 x 10^ points, and we have observed a typical 
statistical uncertainty of about 1% in our results. In Fig. [2] our estimate for 

^1-loop 

is shown as a function of (dimensionless) temperature. Notice the sudden drop to 
zero near Ac = 13.87 which is due to the decoupling of the massive modes aj^^. The 
functional shape is similar to that of the modulus of the leading two-loop correction: 
There is a maximum at A ~ 20 and a very rapid decay to the right of this maximum. 
Notice, however, that the value of the maximum is suppressed by a factor of about 
10^^ as compared to the smallest two-loop correction [10]. In Fig.[3]we present our 
estimate for j^^^^ as a function of (dimensionless) temperature. The maximum of 
this contribution is comparable to the smallest two-loop correction. 

For diagram C we have chosen a compact set {xi, X2, X3 < R, —R < 2/1 < -R, — 1 < 
2^12, ^13 < +1; Z23,i < Z23 < Z23,u} in which the Monte-Carlo method samples points. 
We have varied R in the range 0.1 < i? < 15 and have used samples with up to 
6 X 10® points. We have not found any point which satisfies all of the conditions ( |T7I) . 
This is physically suggestive since the diagram describes annihilation or creation of 
two massive on-shell modes into or out of two massless off-shell modes. Typically, 

^The a;i-integration is performed analytically in order to eliminate the (5-function in the original 
integrand. There are eight zeros of the argument of the (5-function in xi some of which turn out 
to be complex and thus can be discarded. 
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Figure 3: An upper estimate for the modulus of the pressure contribution |APb| 
due to diagram A in Fig.[Tl The plot shows this estimate normalized to the one-loop 
result Pi.ioop, see Eq. (jl]). 

the off-shellness of a massless modes is comparable to the mass of the massive 
mode. This very fact, however, is in stark contradiction with condition ([8]). At loop 
order three we thus have found a first example of an irreducible loop diagram that 
vanishes. We expect that this situation occurs very frequently at higher loop orders, 
see discussion in [9]. 

5 Conclusions 

In this paper we have performed estimates on the moduli of the three irreducible 
three- loop diagrams which contribute to the pressure of a thermalized SU(2) Yang- 
Mills theory being in its deconfining phase. Our results are consistent with the 
general arguments in [9] which imply a rapid convergence of the loop expansion. 
Namely, one of these three diagrams vanishes exactly because the constraints on the 
loop momenta, which emerge in the effective theory, imply that the support for the 
integration is empty. So the situation that an irreducible diagram is precisely zero, 
which in p] is argued to occur at a finite loop order, takes place at three-loops for 
the first time. Moreover, we observe that the modulus of the dominating three-loop 
diagram is comparable to that of the smallest two-loop diagram and by a factor 
of ~ 10~^ suppressed as compared to the dominating two- loop diagram [TU]. The 
modulus of the other nonvanishing, irreducible three-loop diagram is by a factor of 
~ 10^^^ suppressed as compared to the dominating two-loop diagram. 
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